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Abstract 
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We construct a class of systems for which quantum dynamics can be expanded around a mean field 
. approximation with essentially classical content. The modulus of the quantum overlap of mean field states 

naturally introduces a classical distance between classical phase points. Using this fact we analytically 
show that the time rate of change (trc) of two neighbouring classical trajectories is directly proportional 
to the trc of quantum correlations. Coherence loss and nonlocality effects appear as corrections to mean 
field dynamics and we show that they can be given in terms of classical trajectories and generalized 
actions. This result is a first step in the connection between quantum and classically chaotic dynamics 
in the same sense of semiclassical expansions for the density of states. We apply the results to the 
nonintegrable (classically chaotic) version of the N-atom Jaynes-Cummings model. 

; 

I Ever since the conception of Quantum Mechanics the classical limit has been a matter of much debate due 

^ ■ to the profound contrasting differences between the classical and quantal descriptions of the world. Although 

the difhculties in building a bridge between quantum and classical mechanics are well known, we start by 
reviewing the ones which are of relevance for the present contribution. As far as kinematical differences are 
concerned, already at the level of a point particle, striking differences appear. While the definition of the state 
■ of a classical particle is of local character and given by a point in phase space, the quantum counterpart of the 

I definition of a particle state is given by a vector in Hilbert space which cannot simultaneously be ascribed a 

well defined value for position and momentum. The closer one can get to the classical situation arc minimum 
uncertainty wave packets. Quantum states are therefore usually nonlocal. Also, the linear character of the 
Hilbert space has the immediate consequence that superposition of (minimum uncertainty) states are also 
possible states and, in fact, constitute the vast majority of allowed quantum states. The situation gets 
even cloudier when two degrees of freedom are involved: Classically one can always describe a two particle 
state in terms of the coordinates and momenta of each one of them. Quantum mechanically, however, this 
situation only holds if the two particles, initially in a factorized state, do not interact. The Hilbert space 
structure allows for states which cannot be written as a direct product of vectors in the individual Hilbert 
spaces of each degree of freedom. This essentially quantum property is usually named entanglement. Much 
investigation and progress both on the theoretical and experimental sides have been achieved recently, jl] 

From the dynamical point of view one of the essential differences has given rise to an important research 
area nowadays: classical chaos, a phenomenon whose root lies on the nonlinearity of Newton's equation. The 
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relevant question here is how to identify the quantum counterpart of classical chaos. A major step in this 
direction was given by Bohigas and collaborators who conjectured that spectral properties of integrable and 
nonintegrable systems should be very different Q . Thereafter many numerical investigations confirmed such 
conjecture and a few exceptions where found. From the analytical point of view a most relevant formula 
was derived by Gutzwiller connecting level densities of very general quantum systems with classical periodic 
orbits and their actions]^. In what concerns the connections between classical and quantum dynamics it has 
recently been proposed by Zurek that the rate of entropy production can be used as an intrinsically quantum 
test of the chaotic vs. regular nature of the evolution^. Several numerical tests of this conjecture can also 
be found 1^, Anyway, analytic results in this context are scarce. This letter is a step in the direction of 
filling in this gap. 

Given the considerations above we restrict ourselves to the study of the (two degrees of freedom) hermitian 
bilinear hamiltonians 

H = Y^ a.,A. + ^ P,B, + ^ j,^,AB, , (1) 

i 3 i,3 

where i,j = 0,±, A and B are chosen among the generators of either the /i(3) (Heisenberg group) or SU{2) 
groups. Index is associated with the operator a^a {Jz), the index + with the operator (J+) and the 
index — with the operator a (J_) for the algebra h{3) (su(2)). The coefficients a, (3 and 7 are constants or 
given functions of time. 

Although seemingly trivial this class of systems encompasses a rich variety of dynamical behavior in- 
cluding the model we shall use for illustration whose clasical limit is chaotic. The choice of the groups h{3) 
and SU{2) is due to the fact that several of the problems mentioned above can be circumvented at the 
lowest order, the mean field. Due to the bilinear character of the hamiltonian, a time dependent mean field 
solution will be products of coherent states whose labels satisfy the classical limit of Heinsenberg's equations 
of motion as can be easily verified. No quantum dynamical nonlocality effects appear at this level (lowest 
order) unlike mean field approximations for other systems. All quantum corrections will appear in next to 
leading orders as we will show. 



The zeroth order approximation: The mean field approximation (MFA) and the classical limit. 

Our zeroth order approximation is defined in the following way: 

a) We consider here states of the form of products of coherent states and a phase 

|V(0) =exp(^7;(^))|x(^))® |y(t)), (2) 
where the minimum uncertainty coherentes states \z{t)) are defined by \z{t)) = 'D{z{t))\0) with 

-nf m\ / exp[0(t)at-z*(t)a] 

''^^(^^^^l exp[S^^(.(t)J+-z*(t)J_)], (3) 

for h(3) and su(2) respectively, and the fiducial state |0) is the Fock state |n = 0) for h(3) and the Jz 
eigenstate \J,—J) for su(2). 

b) Linear combinations are not allowed as initial conditions (this circumvents problems with the superposition 
principle). Also this MFA for the systems (|l|) will leave invariant the manifold of the chosen set of trial 
functions. 

At this point it is important to mention that the algebra h(3) could be easily enlarge to include of 
harmonic oscillator algebra with little effort, but the extra terms give rise to nonlocality effects (squeezing 
dynamics) already at this lowest order which we would like to avoid. Notice also that generalization to n 
degrees of freedom m-linear hamiltonians (m < n) is straightforward. 

Now we solve the mean field Schrodinger equation 

{zdt-H)\m), (4) 
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where h = 1 and the mean field hamiltonian (MFH) is 

H=1a^ HB{x{t)) + HA{y{t)) ® Ib, 
with HB{x{t)) = {x{t)\H\x{t)) . Using the explicit form of the hamiltonian (|l|) we find 

ilA(t) -^a,(OA, + /^(t), (5) 

i 

where aiit) = Oi + J^jliAvl^jlv) ^nd /^(t) = J2j Pj{y\Bj\y)- The expression for HB{t) is completely 
analogous to the expression for HA{t). Since the terms fA{t), fsit) give rise to global phases we neglect 
them. Due to the structure of the MFH the solutions of eq. (^) will preserve the form (|^), and their time 
dependence will be completely specified by the solutions of the following equations of motion 

i = —iaoz — iaj^ (h(3)) (6) 

i = — m+ — iao2: + j(a+)*z^ (su(2)) (7) 

= {Q\V\z{t)){idt-h{t))V{z{t))\Q). (8) 

Observe that the total phase rj{t) is the sum of the partial phases rjx{t) and Tjy{t). The nonlinearity of these 
equations arise from the self consistency of MFA. If the labels are scaled as z = ^/AJZ for h(3), z = Z for 
su(2), and time as t = ic/(4t/), then the equations for Z(tc) will become independent of J and correspond 
to the classical limit of Heisenberg's equations of motion. From equation (||) it follows that the phases riz{t) 
are generalized actions of the coherent state \z{t)). They are of course absent from the classical limit, but 
will be crucial for the quantum corrections. 

Thus, we have shown that systems with hamiltonian (|l]), in the mean field approximation satisfy all the 
requirements we wanted: labels with classical physical meaning, no superposition principle, minimum quan- 
tum nonlocality effects, and hamiltonian equations for labels which coincide with both the exact Schrodinger 
equation and the Heisenberg equations of motion. 

Quantum corrections. Now we turn to the corrections to the MFA. The exact Schrodinger equation for 
the whole system (with the tilde indicating the Schrodinger picture) 

idt\m) = {H{t)+~m)m)), 

where A(t) = H — H{t), can be written in the MFA interaction picture (MFAIP) as (the absence of the tilde 
indicating the MFAIP) 

zdtim) = m\m) ~ {AmKB^t) - {B.mmt)), (9) 

again up to time dependent terms that give rise to a global phase. Here we have defined A{t) and {A){t) as 

A{t) = U\t,0)AU{t,0), 
{A){t)^{0\vHx{t))AV{x{tm, 

where U{t, 0) is the evolution operator of the MFA, the product of the evolution operators for each degree 
of freedom. The equations for B{t) and {B){t) are completely analogous. Equation (^ possesses a natural 
expansion 

|^t) = (l-i fdtiAt,- fdtJdt2At,At, + ---)\I), (10) 
Jo Ja JQ 
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where |/) is the initial state \x{0),y{0)). Now, let us see that all the terms in this expansion are readily 
written in terms of classical quantities. For example, the first correction term gives 



mt)) = \i)-^c{t)\D) 

where C(t) is a function to be defined below and \D) is a generalized coherent state, orthogonal to the initial 
coherent state, given by \D) = \D{z{0))) = l?(z(0))|l) with the reference state |1) given by the Fock state 
\n = 1) for h(3) and by | J, — J+ 1) for su(2). We immediately see that the first correction already introduces 
all of the effects we avoided at the mean field level: superposition of states, nonlocality and entanglement. 
State \D) is sometimes called a doorway state. In order to get the first order result we made use of the 
following identity which is a result of the present group(s) structure. In fact, it is this relation which enables 
one to express all order corrections in terms of classical trajectories and generalized coherent states. 

AD{x) = D{x) (^gf,{x)Ak + kf{x)^ , (11) 

where gik{x) are functions of x specific to each one of the groups in question. Similar relations hold for the 
degree of freedom B. The coefficient C{t) = c{ti)dti, where 

is clearly only a function of classical trajectories and the corresponding actions, So{t, 0) — dti {'rix{t)+Tiy{t)) 
and Si{t) = Si[x{t)] + Si[y{t)] with 

SMt)] = t dr{l\VHzir)){idr - hir))V{z{r))\l). 



Jo 

It is clear that the remaining corrections can also be written in terms of classical trajectories (and actions), but 
their quantum content will not be as transparent as in the leading correction. In fact, the second correction 
can be written as having a term proportional to the initial coherent state, but also terms proportional to 
generalized coherent states |a;o,yi), \xi,yo), and |x2, 2/2)7 where the subindices refer to the fiducial states. 

Sensitivity to Initial Conditions: A Formal Nonperturbative Result Classically one of the basic 
ingredients to define chaos is the high sensitivity to initial conditions. A formalization of this condition is 
heavily based on the notion of distance between trajectories. In establishing a quantum counterpart of this 
condition it is important to introduce a quantum measure of distance between states. A natural measure is 
given by the square modulus of the scalar product .[| In what concerns our mean field approximation, the 
squared modulus of the scalar product between different states of the manifold allows for a direct association 
of the distance between states with distances between phase space trajectories, since |(2;i|z2)P is given by 

r exp(-|zi-Z2p) forh(3), 

The important quantum tool which allows us to investigate the sensitivity to initial conditions of the 
quantum dynamics and eventually make connection to the well known classical limit is the overlap between 
two time dependent states, which evolved from different initial conditions. It is well known that for unitary 
evolutions the scalar product is conserved in time. Observe that different initial states correspond to different 

'^Some proposals have been made related to scalar product of wavefunction evolved from different hamiltonians (not different 
wavefunctions). See A. Peres, Phys. Rev. A 30, 1610 (1984); R.A. Jalabcrt, and H.M. Patawski, Phys. Rev. Lett. 86, 2410 
(2001). 



4 



MFAs, since the MFA is state dependent, due to self consistency. Thus the scalar product between two 
different initial states is to be written as 



(x'(0),y'(0)|x(0),y(0)) = (x'W,y'WN(t),2;(t)) + (x'W,y'WI('5^cW)^NW,2/W) (13) 
+ {x'{t),y\t)\SQcit)\x{t),yit))+{x\t),y\t)\iS'Qcm^SQc^^^^^ 

where we have written U{t), the exact quantum time evolution operator as U'''\t){l + Sq'(-;), with tj^'^{t) the 

MFA evolution operator corresponding to state \x{tY'\y{tY'^) and 1 + (5q|^ the evolution operator for the 
quantum corrections. 

The first term on the rhs, {x' {t)^y' {t)\x{t)^y{t)) contains the mean field approximation and is given by 

^-d{x{t)-x'(t))/2^t'i>(x(t),x'{t))^-d{y(t)-y'{t))/2^i'i'(y{t),y'{t)) ^ 

where $ is some phase which also depends on the classical trajectory and the corresponding group, and 



functions d can be determined by comparison with (12). This matrix element is proportional to the distance 



between the labels which, in the present case, corresponds precisely to the classical trajectories. Since the 
exact evolution preserves overlap, the sum of this term with the other three of eq. (|l3|), which contain the 
quantum corrections, should be conserved in time. Observe that the rate of change of this overlap can have 
two distinct origins, dephasing and/or change of the modulus. Both changes should be compensated by 
quantum corrections, but only the later one can be unambigously connected to classical chaos, since, if the 
system is classically chaotic this distance will exponentially grow and, as a consequence, the overlap involving 
only the MFA will decrease accordingly. This is a quantum counterpart of the fact that classically chaotic 
systems exhibit high sensitivity to initial conditions. The corresponding quantum system will exhibit a high 
sensitivity to the initial state in what concerns the production rate of non unitary quantum corrections to this 
overlap. We remark that this result is exact and independent of the approximation used. For the argument, 
however it has been crucial to separate the mean field contribution explicitly. It should also be emphasized 
that the time scale for the overlap quantum corrections is essentially linked to the Lyapunov exponents. 
However, other observables will have different time scales, sometimes much shorter, as for example that 
for the entanglement process. Consequently, entanglement is not always a good measure of classical chaotic 
behavior, at least for short times, unless the exponential separation of neighbouring classical trajectories also 
occurs at very short time scales. Should the exponential separation occur at early times, a significant increase 
in linear entropy will be noticed. In the cases the two time scales are very different this effect, although 
present, will be rendered less conspicuous by the time development of quantum correlations stemming from 
the other sources. This will become clear in the example below. 

In order to characterize the degree of entanglement we will calculate the idempotency defect (or linear 
entropy) 5{t) = 1 — TrA(pA(i))^, where the reduced density pA{t) is given by pA{t) = TrB|V'(*))(V'(OI- 
Observe that this measure of entanglement does not depend on the picture used to calculate it. Using the 
expansion ([lo| ) up to second order, writing |/) = |/^) ® \Ib), and calculating the idempotency defect in the 
MFAIP we obtain to second order 

Sit) =AKe [ dti [ dt2 c*{ti) c{t2). 
Jo Jo 

Application to the classicaly chaotic maser model The classically chaotic maser hamiltonian 

G G' 

H ^eJ^+ ijjo^a + —^{a^J^ + aj+) + -^{a^ J+ + aJ-) 
V J y J 

belongs to the class of bilinear hamiltonians (|l]), where the field (atomic) degree of freedom A (B) is related 
to the h(3) (su(2)) algebra. Field coherent states are characterized by x{t) while spin coherent states by 
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y(t). The MFH is of form (||) for each degree of freedom with 

.Gy + G'y 



.Gx + G'x* 



, a_ = (a+)* ao = 
, a_ = (a+)* ao = e, 



for the field and atomic MFH respectively. Also, the equations of motion for each degree of freedom are 
equations (|^) for the field degree of freedom and (0) for the atomic one. For this model it is a simple matter 
to give an analytic expression for c(t) , the key ingredient for evaluation of both the first order correction for 
the state, and the idempotency defect. 

/2pi(So,-j-Si,-j+i)(t) 

where Sq^-^j (S'i._,/=i) is the generalized action for the coherent state with fiducial state \n = 0) (g) \J, ~J) 
{\n = 1) (g) |J, — J + 1)). In this case the doorway state \D) is given by \D) = 

The time development of the magnitude of the overlap in the MFA between two coherent states centered 
in the classically chaotic phase space region and also the overlap between other two states chosen in the 
regular region are shown in fig. and illustrate the effect of the classically chaotic motion on quantum 
dynamics is dramatic. In effect, for the times the overlap changes appreciably, one can check that the classical 
trajectories of the chaotic region in question are also correspondingly well set apart. It is clear that for the 
magnitude of the overlap in MFA the time scale of correlation effects is essentially dictated by classical 
dynamics. Of course, this needs not hold for other quantum observables. Entanglement, for example, is a 
quantum property with a smaller time scale. The expected sudden increase in this quantity at the time the 
modulus of the overlap diminishes, disappears in the midst of contributions of several quantum effects other 
than the one related to the classical limit (see ref. |5j). Our analytical approximation for entanglement breaks 
down for times of the order of 2 for initial conditions both in the chaotic and regular regions. Generalization 
of these results to other quantum systems is possible, but quantum effects such as nonlocality will be already 
present at the lowest order, and other analytical approximations should be advanced. Work along these lines 
is in progress. 
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Figure 1: Squared modulus of the overlap between two neighbouring states \{xi{t),yi{t)\x2{t),y2{t))\'^ 
in the MFA for mean energy E=8.5, J=9/2 in a resonant (e = 1 = w) non integrable case with G=0.5, 
G' = 0.2 for conditions in the chaotic region (continuous lines) and regular region(dashed line). Chaotic 
initial conditions (xi,yi) = (5.7263433,-0.24253563), (^2,^2) = (5.7778567,-0.26845243). Regular initial 
conditions {xi,yi) = (3.615516,0.53452248), {x2,y2) = (3.68977334,0.50086791). 
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